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Intermittent permeation of cylindrical nanopores by water 
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Molecular Dynamics simulations of water molecules in nanometre sized cylindrical channels con- 
necting two reservoirs show that the permeation of water is very sensitive to the channel radius 
and to electric polarization of the embedding material. At threshold, the permeation is intermittent 
on a nanosecond timescale, and strongly enhanced by the presence of an ion inside the channel, 
providing a possible mechanism for gating. Confined water remains surprisingly fluid and bulk-like. 
Its behaviour differs strikingly from that of a reference Lennard- Jones fluid, which tends to contract 
into a highly layered structure inside the channel. 



Highly confined molecular or ionic fluids, including 
fluid films compressed by surface force machines [|l], or 
trapped in nanoporous materials like zeolites 0] or car- 
bon nanotubes || Q, exhibit distinctive structural, dy- 
namical and phase behaviour which may differ widely 
from their corresponding bulk properties. Recent experi- 
ments ^ and simulations [Q demonstrate the interesting 
behaviour of water under conditions of extreme confine- 
ment, which is related to the disruption of the hydrogen- 
bond network between H2O molecules, depending on the 
hydrophobic or hydrophilic nature of the confining sur- 
face. A good understanding of the molecular organiza- 
tion of confined water is essential to the study of the 
selectivity and permeation mechanisms of ion channels 
through membranes Q . 

While much of the published simulation work is on re- 
alistic or semi-realistic models of specific ion channels , 
the present letter examines the generic behaviour of wa- 
ter permeation in a highly simplified model of a channel, 
as a function of a small number of physical parameters 
characterizing the channel. The proper inclusion of end 
effects and Coulombic interactions increases the relevance 
of the present results for real ion channels, compared to 
earlier studies of infinite channels ^, [l(| . 

The model considered in the present work is that of a 
cylindrical channel of radius R and length L through a 
slab of dielectric material, of relative permittivity e (rep- 
resenting the membrane) , which separates two reservoirs 
containing water molecules and eventually ions (repre- 
senting the intra and extra-cellular compartments), sym- 
metrically placed on both sides of the slab (see Fig. [l]). 
In our Molecular Dynamics (MD) simulations JTT|, the 
two reservoirs, which are connected through a common 
x-y plane under periodic boundary conditions, contain 
typically N — 500 water molecules each. The number 
of water molecules in the reservoirs fluctuates somewhat 
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FIG. 1: Schematic representation of a cut through the simu- 
lation cell. The z— axis coincides with the axis of the cylin- 
drical channel connecting the two reservoirs; the interior wall 
region is a continuum of permittivity e. Periodic boundary 
conditions apply in all three directions. The channel radius 
and length are shown at their effective values, R e ff < R and 
L e ff > L, as explained in the text. 



when the cylindrical channel fills or empties; in order to 
maintain a constant density of 0.996 gr cm -3 , the reser- 
voir length parallel to the axis of the channel (taken to be 
z) and the z-components of the molecular position vec- 
tors are scaled regularly using a version of the Berendscn 
barostat Q. 

The water molecules interact via the standard SPC/E 
pair potential [T^ | which involves a Lennard-Jones (LJ) 
potential between the O-atoms (a = 0.3169 nm, e = 
0.6502 kJ moF 1 i.e e/k B = 78.2 K ), and Coulombic in- 
teractions between the sites associated with the O atoms 
(q = -0.8476) and the two H atoms (q = +0.4238). The 
simulations reported in this letter contain no ions or just 
one K + C1 _ pair, which interact with the water molecules 
as in ref . || . The radial interaction between the channel 
surface and the water molecules is taken to be the LJ 
potential between a CH3 group and an O atom |Q , inte- 
grated over the surface of an infinitely long cylinder [fl5| . 
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The same potential, now as a function of the distance 
from the planar wall, is used for the membrane- water in- 
teraction, with an appropriate rounding and shifting pro- 
cedure to avoid any discontinuities. A similar procedure 
is used for ion-surface interactions. 

Since water molecules and ions carry charged sites, 
they polarize the membrane modelled here by the dielec- 
tric slab. The interface between slab, channel or reservoir 
is treated as infinitely sharp (dielectric discontinuity), so 
that the polarization charge is entirely localized on the 
surface of the slab and can be calculated by the vari- 
ational procedure described in ref. ]l6|] ; the numerical 
implementation requires the use of a grid spanning the 
cylindrical surface and the planar interface between the 
slab and the reservoirs. Most of the results presented 
in this letter neglect membrane polarization, i.e. e = 1 
is assumed; the importance of polarization on the be- 
haviour of confined water is assessed later in the letter, 
by comparing data obtained for e = 1 and e = 4. 

The present work is restricted to equilibrium fluctua- 
tions. Starting from initial water configurations with an 
empty channel, the number of molecules within the chan- 
nel is monitored as a function of time, for total times up 
to 6.5 ns. Simulations were carried out for a rather short 
channel of length L = 0.8 nm. The channel radius R was 
varied between 0.4 and 0.7 nm, but in view of the soft na- 
ture of the surface-water interaction, the corresponding 
effective radii of the cylinders accessible to the centres of 
the O atoms (i.e. such that the wall-0 repulsive potential 
is less than ksT) are Reff = 0.17 — 0.47 nm. Similarly 
the effective length of the channel is £ e // — 1-26 nm, 
which is roughly the length of the selectivity filter in the 
KcsA channel, the structure of which has recently been 
resolved by X-ray diffraction [jl7| 

When the channels are filled, the radial density pro- 
files of O and H atoms can be calculated by averaging 
over all water configurations generated in the MD runs. 
In order to gain a better understanding of the impor- 
tance of the hydrogen-bond network within the channel, 
simulations were also carried out for a reference system 
of non- associative molecules obtained by removing the 
H atoms, i.e. an imaginary system of "OW" atoms in- 
teracting solely through the Lennard- Jones part of the 
SPC/E potential. While the H2O simulations were run 
at room temperature, the reference runs for the OW 
atoms were carried out at T = 120 K and a reduced 
density p* — na 3 /V — 0.85, conditions under which 
the bulk LJ system is expected to be liquid. Under 
these conditions we also find liquid-like diffusion in the 
pore, whereas at the lower temperatures T — 56.2 K (i.e. 
T* = k B T/e = 0.72, close to the LJ triple point|l|) and 
T — 85 K the OW atoms form a disordered solid in the 
pore. The number n(t) of OW atoms within channels of 
two different radii are shown in Fig. |as a function of 
time. As expected, there are large fluctuations in n(t), 
but the key observation is that the LJ particles fill the 
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FIG. 2: Number n(t) of OW atoms inside channels of radius 
R = 0.45 nm (R eff = 0.22 nm) and R = 0.6 nm (R eff = 
0.37 nm), in bold, as a function of time (in ps). 

channel for all radii considered. The mean density of 
particles inside the channel, (n(t)) / '(nR^fL) is about 
1.7 times larger than the bulk density in the reservoirs, 
but the confined atoms remain in a fluid state, as char- 
acterized by a diffusion constant along the cylinder axis 
(2 x 10~ 5 cm 2 s _1 ) comparable to the bulk diffusion con- 
stant (6 x 10 -5 cm 2 s _1 ). The radial density profile p(r) 
is highly structured, as shown in Fig. 4d, indicative of 
layering. 

The behaviour changes dramatically when the OW 
atoms are replaced by H2O molecules. The numbers of 
molecules inside the channel are shown as a function of 
time in Fig. ^ for four different cylinder radii. When 
R = 0.55 nm, the channel shows no sign of filling and the 
small fluctuations in n(t) away from zero reflect attempts 
by small numbers of molecules to penetrate the chan- 
nel at either end. At the slightly larger channel radius 
R = 0.6 nm, a strikingly intermittent behaviour is ap- 
parent, reminiscent of recent observations in simulations 
of water in carbon nanotubes |ij and of water in atom- 
istically rough models of ion channels [[|. States where 
the channel is filled alternate stochastically with empty 
channel states, on a time scale of typically Ins. When 
the radius is further increased to R = 0.65 nm, the chan- 
nel appears to remain filled after an initial intermittency, 
although it cannot be ruled out that short periods in the 
empty state could be observed in much longer runs. Fi- 
nally at R = 0.7 nm, the channel appears to remain filled 
throughout, with a mean density of 0.034 A" 3 compara- 
ble to the bulk density in the reservoirs, 0.033 A~ 3 . Dif- 
fusion of the confined molecules along the z— axis occurs 
on a timescale comparable to bulk diffusion (the reser- 
voir diffusion constant and the channel z— axis diffusion 
constant both being 3 x 10 -5 cm 2 s _1 ). Radial density 
profiles of the O and H atoms, averaged over occupied 
states of the channel, are shown in Fig. 4 a-c for three 
values of the radius R. They are surprisingly flat, and 
contrast sharply with the density profile of the Lennard- 
Jones particles plotted in Fig. 4d. The absence of a 
clear-cut layering of the H2O molecules, as well as the 
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FIG. 3: Number n(t) of water molecules inside channels of 
length L — 0.8 nm {L e Jf = 1-26 nm) and radii R — 0.55 nm 
(Reff = 0.32 nm) (a), R = 0.6 nm (R eff = 0.37 nm) (b), 
R = 0.65 nm (R eff = 0.42 nm) (c) and R = 0.7 nm (R eff = 
0.47 nm) (d) as a funtion of time (in ps); note that the total 
durations of the four runs vary; all four runs are for e = 1; 
the inset in frame (a) shows n(t) for the same L and R, but 
with e = 4. 




1 2 3 4 5 

d r/A 



FIG. 4: Radial density profiles p a (r) of O atoms (full curves) 
and H atoms (dotted curves) of water in channels of length 
L = 0.8 nm (L e f f = 1-26 nm) and radii R — 0.6 nm (Reff = 
0.37 nm) (a), R = 0.65 nm (Reff = 0.42 nm) (b) and R = 
0.7 nm (R e ff = 0.47 nm) (c). Averages are taken only over 
times when the channel is filled. Frame (d) shows the radial 
profiles for OW atoms in a channel of the same length and 
R — 0.6 nm (Reff = 0.37 nm), at temperatures T = 85 K 
(dotted curve) and T = 120 K (full curve). 



substantial fluidity of the confined water, characterized 
by the bulk-like diffusion, agree qualitatively with recent 
surface force apparatus measurements of highly confined 
water films ||. 

The results in Fig. || show that the filling (or wetting) 
of a cylindrical channel is very sensitive to its radius, 
with a threshold radius below which water does not pen- 
etrate on this timescale. So far the polarization of the 
confining protein and membrane has been neglected, by 
setting e = 1. One intuitively expects that such polar- 
ization by the charge distribution on the water molecules 
will lower the electrostatic energy, and hence favour the 
wetting of narrow channels. This is confirmed by the 
time trace of n(t) shown in the inset to Fig. ||a from MD 
simulations with e = 4: water molecules are seen to fill 
the R = 0.55 nm channel intermittently, while the same 
channel remains empty in the absence of membrane po- 
larization. The latter effect is thus found to lower the 
threshold radius above which water molecules can pen- 
etrate inside the channel. The effect is expected to be 
significantly enhanced for larger values of e, which may 
be more appropriate for the channel protein and mem- 
brane. 

Returning to the non-polarizable case e = 1, water can 
be made to fill a narrow channel by placing an ion inside, 
fn an MD simulation a K + ion was held fixed on the 
axis in the middle of a channel of radius R = 0.5nm; 
the channel was then found to be rapidly filled by H2O 
molecules with an average occupation number n 6.5. 
Water molecules were found to penetrate channels even 
with R as small as 0.3 nm under these circumstances. 



A close examination of molecular configurations, radial 
number density profiles and mean occupation numbers 
leads to the conclusion that water fills ion channels be- 
yond a critical, e-dependent radius in a surprisingly uni- 
form manner which is intermittent near the threshold. 
As soon as water wets the channel, the water molecules 
appear to fill the available volume , with a mean den- 
sity comparable to the bulk density of liquid water. In 
this series of simulations, no long-lived single- file chains 
of water molecules were observed. Using an energetic 
criterion [^o| the mean number of hydrogen bonds was 
found to be 2.5 in a channel of R = 0.7 nm, compared to 
3.5 in the reservoir. 

In summary we have investigated a highly simplified, 
generic model for the selectivity filter of ion channels, 
which is entirely characterized by three physical parame- 
ters R, L and e, once one has fixed the interaction param- 
eters of the wall. The present investigation was limited to 
L = 0.8 nm. At the threshold radius, water fills the chan- 
nel intermittently, and the fraction of time over which the 
channel is in the occupied state depends sensitively on R 
and e. The timescale for the intermittency is expected 
to increase with channel length L, because larger fluctu- 
ations at the channel ends will be required for channel 
filling. Finite, physiologically relevant concentrations of 
salt will be included in future work. In view of the pre- 
liminary results with one ion fixed within the channel, it 
is expected that permeation of narrow channels by wa- 
ter molecules and ions is a very cooperative effect, where 
the polarization of the channel surface will play a crucial 
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role, which may provide a simple scenario for ion channel 
gating. 
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